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By  using  first-principles  electronic  structure  calculation  and  Boltzmann  transport  equation,  we  investigate 
the  impact  of  gallium  (Ga)  doping  on  the  thermoelectric  property  of  [0001]  zinc  oxide  nanowires 
(Zni_xGaxO  NWs).  Our  results  show  that  the  thermoelectric  performance  of  the  Zni_xGaxO  NWs  is 
strongly  dependent  on  the  Ga  contents.  The  maximum  achievable  room  temperature  thermoelectric  figure 
of  merit  in  Zni_xGaxO  NW  can  be  increased  by  a  factor  of  2.5  at  Ga  content  of  0.04,  compared  with  the 
ZT  of  pure  ZnO  NWs.  This  may  open  up  ZnO  NW  arrays  applications  in  possible  thermoelectric  energy 
generator  and  cooler. 

©  2011  Elsevier  B.V.  All  rights  reserved. 


Recently,  zinc  oxide  (ZnO)  nanowires  (NWs)  have  shown  promise 
for  nanodevice  applications.  Through  the  coupling  of  piezoelectric 
and  semiconducting  properties  of  ZnO,  nanoscale  mechanical  en¬ 
ergy  can  be  converted  into  electrical  energy  by  using  ZnO  NW 
arrays.  Moreover,  compared  to  traditionally  Ti02  films,  ZnO  NW 
arrays  are  considered  promising  anode  material  for  photovoltaic 
device  because  of  the  highly  controllable  single-crystalline  mor¬ 
phology  which  provides  direct  electron  transport  pathway.  In  ad¬ 
dition  to  piezoelectric  [1]  and  optoelectronic  applications  [2-4], 
thermoelectric  effect  is  also  an  important  approach  to  the  solu¬ 
tion  to  the  energy  crisis  by  converting  waste  heat  into  electricity. 
The  energy  conversion  efficiency  of  a  thermoelectric  material  is  de¬ 
scribed  by  the  dimensionless  figure  of  merit  ZT,  which  is  a  function 
of  Seebeck  coefficient,  electrical  conductivity,  absolute  temperature 
and  thermal  conductivity.  The  study  of  the  thermoelectric  prop¬ 
erty  of  ZnO  ceramics  has  been  motivated  by  the  reports  of  the 
improved  electrical  conductivity  and  reduced  thermal  conductivity 
by  employing  elements  (Al,  Ga,  Mn)  as  dopants  [5-7].  In  the  A1  and 
Ga  dually  doped  ZnO  ceramics,  ZT  is  about  twice  of  that  of  pure 
ZnO,  which  can  arrive  0.65  at  1247  K  [5].  A  natural  question  comes 
promptly:  How  to  improve  the  thermoelectric  property  of  ZnO 
NWs?  Although  some  research  works  have  been  done  on  the  elec¬ 
trical  properties  of  ZnO  NWs  [8-11],  very  little  attention  has  been 
paid  to  the  thermoelectric  properties  of  ZnO  NWs  [12].  In  this  Let¬ 
ter,  we  will  study  thermoelectric  properties  through  first-principle 
calculation,  ID  Boltzmann  transport  theory  (BTE),  and  molecular 
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dynamics  simulations.  In  particular,  we  will  explore  the  impacts  of 
Ga  doping  on  the  thermoelectric  performance  of  Zni_xGaxO  NWs. 

We  focus  on  Zni_xGaxO  NWs  oriented  along  the_[0001]  direc¬ 
tion,  with  hexagonal  cross  section  enclosed  by  (1010)  facets,  be¬ 
cause  this  is  the  commonest  used  surface  structure  in  experiments 
and  simulations  [8-10].  It  has  been  demonstrated  that  with  the 
NW  diameter  increases,  the  thermoelectric  power  factor  decreases, 
while  thermal  conductivity  increases,  which  results  in  the  reduc¬ 
tion  in  thermoelectric  figure  of  merit  [13].  Moreover,  larger  system 
size  is  challenge  to  current  computational  capability.  As  the  main 
concern  of  this  work  is  the  impact  of  Ga  doping  on  the  thermo¬ 
electric  property  of  Zni_xGaxO  NWs,  only  thin  ZnO  NWs  with  the 
diameter  of  about  0.7  nm  are  considered  here.  The  atomic  struc¬ 
ture  of  Zni_xGaxO  NW  is  initially  constructed  from  bulk  wurtzite 
structure  of  ZnO  solid.  The  Zn  atoms  are  randomly  substituted  by 
Ga  atoms  and  the  geometry  is  relaxed  to  its  closest  minimum  to¬ 
tal  energy.  A  supercell  approach  is  adopted  where  each  wire  is 
periodically  repeated  along  the  growth  direction.  The  size  of  the 
supercells  in  the  transverse  plane  is  large  enough  (>  15  A  from 
surface  to  surface).  Two  repeated  units  along  the  growth  direction 
are  required  for  Zni_xGaxO  NWs  and  the  relaxed  supercell  length 
along  the  wire  axis  is  about  10.4  A  [14].  The  total  atom  number 
of  pure  ZnO  NWs  is  52,  with  26  Zn  atoms  and  26  O  atoms.  In  this 
Letter,  the  density  functional  theory  (DFT)  calculations  are  carried 
out  by  DMol3  package  [15].  The  DFT  calculations  are  performed  by 
using  generalized  gradient  approximation  with  the  functional  pa¬ 
rameterized  by  Perdew,  Burke,  and  Ernzerhof  [16].  And  the  double- 
numerical-polarization  basis  set  is  employed.  Self-consistent  field 
calculations  are  done  with  a  convergence  criterion  of  10-6  hartree 
on  the  total  energy.  The  Brillouin  zone  integration  is  performed  us¬ 
ing  a  1  x  1  x  1 6  Monkhorst-Pack  k-point  grid.  All  the  structures  are 
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Fig.  1.  The  geometry  of  relaxed  Zni_xGaxO  NWs  for  x  —  0.04.  (a)  Top  view  and  (b) 
side  view.  Ga,  Zn  and  0  atoms  are  represented  in  blue,  gray  and  red,  respectively. 
(For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  Letter.) 

fully  optimized  with  a  convergence  criterion  of  0.002  hartree/A  for 
the  forces  and  0.005  A  for  the  displacement.  A  real-space  cut-off  of 
4.0  A  for  the  atom-centered  basis  set  is  chosen  to  increase  compu¬ 
tational  efficiency  while  not  significantly  affecting  the  magnitude 
of  inter-atomic  forces  or  the  total  energies.  The  Gaussian  smearing 
of  electron  density  is  applied  with  the  energy  range  of  0.1  eV. 

The  electrical  conductivity  a ,  the  thermal  conductivity  due  to 
electrons  ke ,  and  the  Seebeck  coefficient  5,  are  obtained  from  the 
electronic  structure  with  the  solution  of  one-dimensional  Boltz¬ 
mann  transport  equation  as  described  in  Refs.  [13,17].  The  electron 
relaxation  time  is  an  important  factor  for  the  electrical  transport. 
In  doped  ZnO  thin  film,  the  scattering  process  is  dominated  by 
ionized  dopant  atoms  and  boundary.  An  analytical  expression  for 
the  mobility  of  doped  ZnO  thin  film,  taking  into  account  the  non- 
parabolicity  of  the  band  edge,  has  been  given  by  Ellmer  and  Mien- 
tus  [18].  In  doped  ZnO  nanowire,  there  is  additional  scattering 
mechanism  from  nanowire  surface.  However,  as  experimental  mo¬ 
bility  data  for  single  crystalline  ZnO  nanowires  are  not  available, 
here  we  use  the  experimental  data  of  Zni_xGaxO  thin  film  [18,19] 
to  calculate  the  carrier  concentration  dependent  relaxation  time, 
based  on  the  assumption  of  a  similar  boundary  scattering  mech¬ 
anism  in  thin  film  and  nanowire.  Only  perfect  NW  is  considered 
here,  so  there  is  no  defect  scattering.  The  carrier  concentration  is 
defined  as:  n  =  D(E  —  /x)/(F  -  jlc)  dE,  where  f(E )  is  the  Fermi 
distribution  function. 

Fig.  1  shows  the  atomic  structure  for  the  Zni_xGaxO  NW  for 
x  =  0.04  with  the  replacement  of  one  Zn  atom  in  the  supercell. 
x  =  0.08  and  0.12  corresponds  to  two  and  three  Zn  atoms  are  sub¬ 
stituted  with  Ga  atoms,  respectively.  Fig.  2  shows  the  total  density 
of  states  (DOS)  of  Zni_xGaxO  NWs.  For  pure  ZnO  NW  ( x  =  0),  the 
O  2s  states  are  located  around  -18  eV,  which  have  little  inter¬ 
action  with  the  other  valence  band.  The  valence  band  consists  of 
two  main  states,  Zn  3d  and  O  2p  orbitals.  The  Zn  4  s  states  dom¬ 
inate  the  conductance  band.  With  the  Ga  content  (x)  increase,  on 
the  one  hand,  there  is  charge  transfer  from  Zn  4  s  states  to  O  2p 
states,  thus  results  in  the  shift  of  the  center  of  gravity  of  the  DOS 
at  the  O  sites  towards  the  low  energy  regions  and  the  Fermi  en¬ 
ergy  level  into  the  conduction  band  (see  Fig.  2(c)),  which  indicates 
the  n-type  characteristics  of  Zni_xGaxO  wires.  And  on  the  other 
hand,  the  dopant  states  (Ga  4  s  states  and  Ga  3d  states),  appear 
following  the  0  2s  states.  Similar  phenomena  can  be  observed  in 
Ga-doped  bulk  ZnO  [20]. 

Fig.  3(a)  shows  the  Ga  content  (x)  dependent  electronic  band 
gap  and  carrier  concentration.  It  is  obvious  that  for  each  Ga  con¬ 
tent,  the  band  gap  values  depend  on  the  detailed  atomic  structures. 
In  order  to  reduce  the  fluctuation,  the  results  are  averaged  by  over 
five  realizations  such  as  the  case  in  Sii_xGex  NWs  [17].  With  Ga 
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Fig.  2.  The  total  DOS  for  Zni_xGaxO  NWs  for  (a)  x  =  0;  (b)  0.04  and  (c)  0.08.  The 
Fermi  energy  is  set  to  0.  The  dashed  magenta  line  is  used  to  guide  the  eyes.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred 
to  the  web  version  of  this  Letter.) 


content  (x)  increase,  the  averaged  band  gap  decreases.  This  is  con¬ 
sistent  with  the  experimental  observed  energy  band  gap  in  [21].  In 
the  following  calculation  of  the  thermoelectric  properties,  the  DOS 
of  specific  NW  structure  is  used,  which  is  most  stable  with  the 
lowest  total  energy.  The  carrier  concentration  increases  as  the  Ga 
content  increase.  With  x  =  0.08,  the  carrier  concentration  is  about 
6.0  x  1020  cm-3,  and  n  =  2.0  x  1021  cm-3  when  x  =  0.12.  How¬ 
ever,  the  maximum  electron  concentration  can  be  achieved  in  ZnO 
NWs  via  doping  method  is  in  the  order  of  1020  cm-3  [22],  thus  the 
maximum  Ga  doped  content  (x)  we  consider  in  this  Letter  is  0.08. 
Fig.  3(b)  shows  the  Ga  content  dependence  of  a  and  S.  a  increases 
and  S  decreases  as  the  Ga  content  increases  because  there  is  more 
charge  concentration.  At  low  carrier  concentration  range,  the  See¬ 
beck  coefficient  is  about  400  pV/K,  which  is  in  good  agreement 
with  the  experimental  observation  in  Ref.  [12].  Fig.  3(c)  shows  the 
power  factor  (P  =  S2a)  versus  the  Ga  content.  There  is  an  optimal 
Ga  content  (x  =  0.04,  n  —  1.0  x  1017  cm-3)  yielding  the  maximum 
attainable  value  of  power  factor  (Pmax  =  908  pW/mK2). 

ZT  is  one  of  the  central  physical  quantities  that  determine 
the  efficiency  of  thermoelectric  application.  In  the  calculation 
of  ZT,  both  electrons  and  phonons  contribute  to  the  total  ther¬ 
mal  conductivity.  To  calculate  the  phonon  thermal  conductivity 
of  Zni_xGaxO  NW  by  using  molecular  dynamics  simulations,  the 
potential  that  describes  the  interactions  between  the  atoms  is 
needed.  However,  there  is  no  consensus  in  the  literature  about 
the  potential  between  Ga,  O  and  Zn  atoms.  In  the  Zni_xGaxO  NW, 
compare  to  the  original  Zn  atom,  the  doped  Ga  atom  has  differ¬ 
ent  atomic  mass  and  different  interaction  with  the  neighboring 
atoms.  Through  the  investigation  of  thermal  conductivity  by  us¬ 
ing  molecular  dynamics  simulation,  it  has  been  demonstrated  that 
the  difference  in  atomic  mass  has  larger  impact  on  thermal  con¬ 
ductivity  than  the  different  interaction  does  [23].  Thus  here  we 
use  the  Buckingham  potential  [24]  to  describe  the  atomic  interac¬ 
tion  in  Zni_xGaxO  NW,  which  includes  anharmonicity  to  describe 
the  phonon-phonon  scattering  process.  In  the  molecular  dynam- 
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Fig.  3.  (a)  The  averaged  electronic  band  gap  and  carrier  concentration  versus  Ga  contents  for  Zni_xGaxO  NWs.  (b)  Electrical  conductivity  o  and  Seebeck  coefficient  S  versus 
Ga  content  x.  (c)  Power  factor  of  Zni_xGaxO  NWs  versus  Ga  contents,  (d)  Phonon  thermal  conductivity  versus  doping  contents  at  room  temperature.  Here  kq  is  the  thermal 
conductivity  of  pure  ZnO  NW. 


ics  simulation,  there  are  20  unit  cells  in  the  longitudinal  direction 
which  corresponds  to  a  length  of  10.4  nm,  and  the  wire  diameter  is 
0.7  nm,  the  same  as  that  in  the  DFT  calculation.  In  order  to  estab¬ 
lish  a  temperature  gradient  along  the  NW  longitudinal  direction, 
the  two  ends  of  NWs  are  put  into  Nose-Hoover  heat  bathes  with 
different  temperature.  The  simulations  are  performed  long  enough 
to  allow  the  system  to  reach  a  non-equilibrium  steady  state.  All  re¬ 
sults  given  in  this  Letter  are  obtained  by  averaging  about  1  x  107 
time  steps,  a  time  step  is  set  as  0.8  fs.  Free  boundary  condition  is 
used  to  atoms  on  the  outer  surface  of  the  NWs.  The  thermal  con¬ 
ductivity  is  calculated  from  ic  =  —  Jl/VT,  where  the  heat  current 
Jl  along  the  longitudinal  direction  is  defined  as  the  energy  trans¬ 
ported  along  the  NW  in  unit  time  through  the  unit  cross-section 
area,  and  VI  is  the  temperature  gradient. 

Fig.  3(d)  shows  that  the  phonon  thermal  conductivity  decreases 
as  the  Ga  content  increases.  Compared  to  the  60%  thermal  conduc¬ 
tivity  reduction  (with  x  =  0.1)  in  Sii_xGex  NW,  the  thermal  con¬ 
ductivity  reduction  in  Zni_xGaxO  NW  is  only  about  10%,  due  to  the 
small  difference  in  atomic  masses  between  Zn  and  Ga  atoms  [25]. 

Combined  the  calculations  of  the  power  factor  and  the  phonon 
thermal  conductivity,  we  show  the  Ga  content  dependent  ZT  in 
Fig.  4.  The  value  of  ZT  firstly  increases  with  the  Ga  content,  reaches 
a  maximum  value  at  x  =  0.04,  which  is  2.5  times  larger  than  that 
of  pure  ZnO  wires,  and  then  decreases.  This  phenomenon  can  be 
explained  as  follows.  ZT  is  contributed  by  both  power  factor  and 
thermal  conductivity.  With  Ga  content  increase,  the  power  fac¬ 
tor  increases,  reaches  a  maximum  value  at  x  =  0.04  and  then 
decreases.  The  curve  of  the  thermal  conductivity  decreases  with 
the  Ga  content  increases.  Thus  when  0  <  x  <  0.04,  ZT  increases 
with  the  Ga  content.  However,  when  0.04  <  x  <  0.08,  the  reduc¬ 
tion  in  the  thermal  conductivity  cannot  offset  the  reduction  in  the 
power  factor,  thus  results  in  the  decrease  in  ZT  with  the  Ga  con¬ 
tent.  It  is  worth  to  mention  that  the  optimal  carrier  concentration 
yielding  the  maximum  attainable  value  of  ZT  is  1.0  x  1017  cm-3. 
This  low  carrier  concentration  suggests  that  the  high  thermoelec¬ 
tric  performance  of  Ga-doped  ZnO  NW  can  be  realized  by  using 
reliable  fabrication  technology.  It  is  worth  pointing  that  the  op- 


Fig.4.  ZTZni_xGaxo/ZTzno  versus  the  Ga  content  for  Zni_xGaxO  NW.  The  solid  line  is 
used  to  guide  the  eyes. 

timal  Ga  content,  4%  yielding  a  2.5  times  enhancement  in  ZT  is 
an  approximated  value  due  to  the  limited  supercell  size.  To  get  a 
more  accurate  optimal  Ga  content,  larger  supercell  size  is  required, 
which  is  challenge  to  current  computational  capability.  However, 
our  results  demonstrate  that  at  least  a  2.5  times  enhancement 
in  ZT  can  be  achieved  under  a  reliable  doping  concentration.  As 
the  Seebeck  coefficient  is  similar  to  the  value  of  pure  ZnO  NW, 
it  is  obvious  that  the  enhancement  in  ZT  is  mainly  due  to  the 
increase  in  electrical  conductivity  and  reduction  in  thermal  con¬ 
ductivity. 

In  summary,  we  have  investigated  the  Ga  content  dependent 
thermoelectric  properties  in  Zni_xGaxO  nanowires.  It  is  found  that 
the  thermoelectric  performance  is  strongly  dependent  on  the  Ga 
content.  The  Ga  doping  has  reduced  the  thermal  conductivity  of 
ZnO  NW,  down  to  less  than  90%  at  Ga  content  of  0.08.  The  max¬ 
imum  attainable  ZT  corresponds  to  an  optimal  Ga  content  of  x  = 
0.04,  which  is  2.5  times  of  that  of  the  pure  ZnO  wires.  Our  work 
open  up  possible  ZnO  NW  arrays  applications  in  thermoelectric  en¬ 
ergy  generation  and  refrigeration. 
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